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Abstract 
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1. INTRODUCTION 



Recently, the observations of high redshift type la supernovae, the surveys of clusters 
of galaxies [I] [I], Sloan digital sky survey (SDSS) [5J and Chandra X-ray observatory jH] 
reveal the universe accelerating expansion and that the density of matter is very much less 
than the critical density. Also the observations of Cosmic Microwave Background (CMB) 
anisotropies |7| indicate that the universe is flat and the total energy density is very close 
to the critical one [8]. The observations though determines basic cosmological parameters 
with high precisions and strongly indicates that the universe presently is dominated by a 
smoothly distributed and slowly varying dark energy (DE) component, but at the same 
time they poses a serious problem about the origin of DE |9J. A dynamical equation of 
state ( EoS) parameter that is connected directly to the evolution of the energy density in 
the universe and indirectly to the expansion of the Universe can be regarded as a suitable 
parameter to explain the acceleration and the origin of DE jTU] [H] - In scalar-tensor theories 
[T2] [HI], interaction of the scalar field with matter ( for example in chameleon cosmology) 
[20]~[25J or with geometry in Brans-Dicke (BD) cosmological models [2B] [30] can be used 
to interpret the late time acceleration. 

On the other hand, by using the well-known geometric variables, Hubble parameter and 
deceleration parameter together with the new geometrical variables, the cosmological diag- 
nostic pair {r, s} ( or statefinder parameters) |31] . the acceleration expansion of the universe 
and differentiation among the cosmological models can be explained in order to better fit 
the observational data. The importance of the statefinder parameters to distinct DE cos- 
mological models is best realized, in particular, when considering the increased accuracy of 
the observational data during the last few years and generality of the DE models. These 
parameters, in a natural next step beyond the well known geometric variables, are to dif- 
ferentiate the expansion dynamics with higher derivatives of the scale factor and to explore 
a series of DE cosmological models, including A cold dark matter (LCDM), quintessence, 
coupled quintessence, Chaplygin gas, holographic dark energy models, braneworld models, 
and so on [32]-[35]. Moreover, since the cosmic acceleration affects the expansion history 
of the universe, to understand the true nature of the driving force, mapping of the cosmic 
expansion of the universe is very crucial [36]. Hence, one requires various observational 
probes in different redshift ranges to understand the expansion history of the universe. One 
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of these tests is the difference in distance modulus measurement of type la supernovae that 
helps us to testify the cosmological models. 

In this paper, we investigate the dynamics of the scalar-tensor theory with tachyon po- 
tential where an scalar field function F(<p) coupled to the curvature and matter lagrangian. 
The model allows scalar field function that might be light on cosmological scales, to couple 
to matter much more strongly than gravity does, and still satisfies the current experimental 
and observational constraints. The cosmological value of such a field evolves over Hubble 
time-scales and could potentially cause the late-time acceleration of our universe. 

The approach we used in this paper is simultaneously solving the field equations and 
best fitting the model parameters and initial conditions with the most recent observational 
date for distance modulus using chi-squared method. Sec. two is devoted to a detailed 
formulation of the cosmological model. In sec. three we solve the system of field equations 
by simultaneously best fitting the model parameters and initial conditions with the obser- 
vational data for distance modulus. In Sec. four, we examine the behavior of the best fitted 
effective EoS parameter of the model and also perform a statefinder diagnostic for the model 
and analyze the evolving trajectories of the model in the statefinder parameter plane. In 
Sec. five, we present summary and remarks. 



2. THE MODEL 



The model is presented by the action, 

S = j [F{^)R - V(<p)^l - <p^# + F^)L m ]^-gd A x (1) 

where R is the Ricci scalar, G is the Newtonian constant gravity, and the second term in the 
action is the tachyon potential. The F((p) and V(ip) are analytic function of the scalar field. 
The last term in the Lagrangian brings about the nonminimal interaction between the matter 
and the scalar field. The variation of action ([!]) with respect to the metric tensor components 
in a spatially flat Friedmann-Robertson- Walker cosmology yields the field equations: 

F[2H + 3H 2 } + 2[p^]H + + ^] - ±Vy/l^fi + ip m = 0, (2) 

6FH 2 + 6[<pf-]H + \{v 2 ) - -j¥= - Pm = (3) 
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where e = 1 — 37 and 7 is the EoS parameter for the matter filled the universe. Variation 
with respect to the scalar field gives 

£ + (3^ + — )(1 - ^ 2 ) = +-^(1 - </? 2 ) 5 + ^(1 - v 2 )> (4) 

where prime is derivative with respect to the scalar field. From the above field equations, 
conservation equation is obtained as 

F 

p' m + 3H(1 +7)p OT = -Pm-p (1 + e) (5) 



which readily integrates to yield, 



M 

P™ ~ a 3(l+ 7 ) F (l+e) ( 6 ) 



with M is a integration constant and eq ^ is the energy constraint corresponding to the 
(0,0)-Einstein equation. In the following we assume that the matter presented in the universe 
is cold dark matter, 7 = 0. We now study the structure of the dynamical system by 
introducing the following dimensionless variables, 

X = ^ Y = t, Z^H, U = fe (7) 

In the following we consider two cases: 

Exponential functions: F(4>) oc e a</) and V(<f>) oc eP^ 

Using equations (^-Q, the evolution equations of these variables become, 

x , _ V-gXY _ 2X | (g) 

Y' = -3Y(1 - V) - W z y2) + ^(6^ + 12)(1 - Y*)i + fa(l - Y*)i^ (9) 
z' = z§- 2 (10) 

^a-^-^l (11) 

where prime in here and from now on is taken to be derivative with respect to N = ln(a) 
and 

# 2 (XZ 2 + 3a 2 (l-F 2 )i) L 2 2 1 ; 2 1 ; 

- 3ct 2 £/(l - Y 2 Y ~ 6" (1 - ^ Y 2 + — ] 



We also have the constraint equation in terms of the new variables as 



1 + ^r t~~~ U = (13) 



By imposing the constraint equation, (13), equations fl8h-(ll) reduce to 



X' = (/3 - a)XYZ - 2X^- (14) 

ri z 

Y = -3Y(1 - F 2 ) - + ^(Q§- 2 + 12)(1 - r 2 )l (15) 

t 1 + sY_ x ^ 



+ 6a(i-y 2 )i— z 



xz 

Z ' = Z §~ 2 (16) 
Power law functions: F(<p) oc </> a and V(</>) oc (p 13 

By taking the scalar function F(<j>) and potential V(4>) in power law form the dynamical 
equations become 



JC= ^- a)X r-2X* (17) 

z(xz 2 )^ h 2 1 ; 

Y' = -3Y(1 - F 2 ) - ^ -f) + Q (6-^ + 12)(1 - r 2 )§ (18) 
+ 6a(l-r 2 )i- ( 



Z(XZ 2 )^ H 2 



Z' = ^4 (19) 
rf = - W -J?*°- W * (20) 



where 



# 2 [X^(X^) + 3a2(l-y2)i] L 2 1 ; 2 1 

+ ^^"^ - 3a 2 f/(l - F 2 )§ - 6« 2 (1 - F 2 )f 

a(a-l)YX X 2 Z 2 (1 -F 2 )i(X 2 zA) 
2 + 4 J 

The constraint Friedmann equation now becomes 



Using the constraint (21), the above equations now reduce to the following equations: 



(g - <*)*Y 9Y H 

i ZJi ~Ejv 

Z(XZ 2 )J^ H 2 

m-Y 2 ) a H 

— H — 6 777 

Z{XZ 2 )J^ XZ(XZ 2 )^ H 



X = 1„„1^ (22) 
y' = -3y(l - Y 2 ) - fJyl ~ ± I + -(6-^ + 12)(1 - y 2 )i (23) 



1 



aY X 



+ Qa(l-Y 2 )l Z{XZ ^ ^ 

XZ(XZ 2 )T^ 

Z' = Z^L (24 ) 

In the following we solve the dynamical systems on both cases by best-fitting the model 
parameters and also initial conditions with the observational data. 



3. OBSERVATIONAL BEST FIT BY USING DISTANCE MODULUS, n(z) 

The difference between the absolute and apparent luminosity of a distance object is given 
by, n(z) = 25 + 51og 10 di(z) where the Luminosity distance quantity, (1l(z) is given by 

W ,) _(! + ,) (25) 

With numerical computation, we solve the system of dynamical equations for X and Y 
and Z in both in both power law and exponential cases. In addition, to best fit the model 
parameters and initial conditions we use the following three auxiliary equations for the 
luminosity distance and the hubble parameter 

dH H . , 

In = H ^ < 26 > 

= (27) 
dN H K J 



z 



-l + e~ N . (28) 



From numerical computation one can obtain H(z) which can be used to evaluate fi(z). To 
best fit the model for the parameter a and (3, the initial conditions Y(0), X(0), Z(0) with 
the most recent observational data, the Type la supernovae (SNe la), we employ the x 2 
method. We constrain the parameters including the initial conditions by minimizing the \ 2 
function given as 

xiN M, m Y( 0) , zm = £ m^M^mm^l , (29) 



where the sum is over the SNe la sample. In relation (j29h , /i* he and /i° 6s are the distance 



modulus parameters obtained from our model and from observation, respectively, and a 
is the estimated error of the n° bs . From numerical computation, Table I shows the best 
best-fitted model parameters in both cases. 

TABLE I: Best-fitted model parameters and initial conditions. 



parameters 


Ct 


P 


X(0) T(0) Z(0) 


Amm 


Exponential 


— 0.23 


1.17 


5.5 -0.6 0.91 


544.530862 


Power law 


— 2.97 


-1.82 


1 0.96 0.91 


546.1139533 



Figures 1-4 shows the constraints on the parameters a and (5 at the 68.3%, 95.4% and 
99.7% confidence levels in both cases of exponential and power law functions. 




a p 

Fig. 1: The graph of the one dimensional likelihood distribution for 
parameters a and /3 in exponential case 




Fig. 2: The graph of the two dimensional likelihood distribution and confidence level for 

parameters a and (3 in exponential case 
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Fig. 3: The graph of the one dimensional likelihood distribution for 
parameters a and f3 in power law case 




a 




-1.92 -1.82 -1.72 -1.62 -1.52 



P 

Fig. 4: The graph of the two dimensional likelihood distribution and confidence level for 

parameters a and (3 in power law case 

The distance modulus, fi(z), plotted in Fig. 5, in both CclSGS, £1X6 best fitted with the 
most recent SNe la observational data for the model parameters and initial conditions using 
X 2 method. 
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z z 

Fig. 5: The best-fitted distance modulus /jl(z) plotted as function of redshift for F(<p) and V(4>) 

left) in exponential form, right) in power law form 



4. COSMOLOGICAL PARAMETERS 

In order to understand the behavior of the universe and its dynamics we need to study the 
cosmological parameters. We have best fitted our model with the current observational data 
by the distance modulus test. The cosmological parameters analytically and/or numerically 
have been investigated by many authors for variety of cosmological models. Simultaneously 
best fitting the model with the observational data gives us a better understanding of the 
solutions and the dynamics of these parameters. Among cosmological parameters, the effec- 
tive EoS parameter and statefinders are given by u> e // = — 1 — 2H/3H 2 , r = H/H 3 — 3q — 2 
and s = (r — l)/3(g — |) discussed here, where q in r and s is the deceleration parameter 
and in r in terms of new dynamical variables for exponential and power law cases can 
be obtained by taking derivative of H. 

In Fig. 6, the effective EoS parameters are shown in both cases of exponential functions 
(left panel) and power law functions (right panel). In both cases the effective EoS parame- 
ters are best fitted for the model parameters and initial conditions with the observational 
data. From the graph, for the exponential case, the EoS parameter approaches zero in 
high redshift as expected for matter dominated universe. It also shows that at about 
z < 0.5 in the past where u e ff < —1/3 the universe begins to accelerate which is consistent 
with observational evidence. The result does not show phantom crossing in the past 
and present. In the case of power law functions the result is very disappointing. The 
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effective EoS parameter become tangent to the phantom divide line in the past which is 
not observationally justified. Note that in both cases the trajectories are best fitted for the 
model parameters and initial conditions with the observational data. The satisfactory in 
the exponential case and discrepancy in power law case demonstrate an advantage for the 
exponential behavior of these functions. 
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Fig. 6: The best-fitted effective equation of state parameter, Lo e ffi plotted as function of 
redshift for F((p) and V(<fi) left) in exponential form, right) in power law form 

Fig 7 shows the best-fitted trajectories of the statefinder diagrams {r, s} in both expo- 
nential and power law cases. From the graph it can be seen that both best-fitted trajectories 
passed LCDM state with {r, s} = {1,0} sometimes in the past. The current value of the 
best fitted trajectory in exponential case and its location with respect to the LCDM state 
can also be observed in the {r, s} diagram. However, in power law case, since the current 
location of the trajectory is far from LCDM state, in order to find its location and value 
in Fig 8 we depict the corresponding dynamical behavior of the satefinder {r, s} against 
N = — ln(l + z) in both cases of exponential and power law. From Fig 7, we see that the 
current value of the statefinder {r, s} in exponential case is {0.48, 0.19}, close to the LCDM 
state. From Fig 8 in power law case it has been shown that the current value of trajectory is 
{—42, —25}, relatively far from LCDM in statefinder diagram. This analysis shows another 
advantage of considering exponential behavior for the functions F(<fi) and V(4>) over power 
law behavior. This is because the recent observational data confirms LCDM as the standard 
cosmological model in the current epoch. 
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Fig. 7: The best-fitted statefinder parameter {r, s} for F(<f) and V(4>) 
left) in exponential form, right) in power law form 
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Fig. 8: The dynamics of the statefinder parameters r and s with respect to N = —ln(l + z) for 
F((p) and left) in exponential form, right) in power law form 

We also reconstructed the best-fitted function in both exponential and power law cases. 
From Fig 9, we see that the trajectory for the best fitted model parameters in exponential 
case shows monotonic increasing behavior in redshift, whereas in power law case shows a 
constant behavior. Both cases have similar behavior in the near past up to the redshift 
z ~ 0.9. 
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Fig. 9: The best-fitted reconstructed ^plotted as function of redshift for F(<f) and V(<p) 
left) in exponential form, right) in power law form 

5. SUMMARY AND REMARKS 

This paper is designed to study the dynamics of scalar tensor cosmology with tachyonic 
potential and non minimally coupled scalar field function with the curvature and matter 
lagrangian. The model characterized by the scalar field functions F((j>) and V{4>). We 
investigate two forms of the exponential and power law behavior for these functions. To 
solve the field equations, in a new approach, we simultaneously best-fit the model parameters 
and initial conditions with the observational data using \ 2 method. The advantage of this 
approach is that the obtained solutions are observationally verified and thus physically more 
promising. 

We then study the cosmological parameters such as effective EoS and statefinder param- 
eters for the model in terms of the best-fitted model parameters and initial conditions. The 
result shows that the best-fitted effective EoS parameter in case of exponential function, 
while does not cross the cosmological divide line in the near and far past, exhibits an obser- 
vationally verified behavior in the past (oJ e ff — 0, matter dominated universe) . In addition 
its current value is within the range of observationally accepted values. On the contrary, 
in power law case, the dynamics of the best-fitted parameter is very unsatisfactory as the 
result shows an accelerating universe at higher redshifts. The best fitted statefinder param- 
eters show that in exponential case the current state of the universe is very close to LCDM, 
whereas in the power law case it is relativity far from LCDM in the statefinder diagram. As 
noted, this behavior shows another advantage of exponential behavior of the functions over 
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power law one. Finally, the best fitted scalar field function, 0, in exponential case displays 
an interesting behavior, smoothly decreasing function 
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